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Introduction 

The Stieltjes (or generalized Euler) constants [2, 6, 10, 16, 20, 21, 23, 27, 29, 30] 
7fc(a) appear as expansion coefficients in the Laurent series about s = 1 of the Hurwitz 
zeta function ((s, a), a generalization of the Riemann zeta function ((s) [3, 11, 16, 18, 
24, 28]. We present both new asymptotic and exact expressions for these and other 
fundamental mathematical constants. 

In the following, we let s(n,m) denote the Stirling numbers of the first kind and 
Cfc(a) = 7fc(a) — (\n k a) /a. By convention, 7^ represents 7fc(l). With B n (x) the 
Bernoulli polynomials, their periodic extension is denoted P n (x) = B n (x — [x]). As is 
customary, Bj denotes the Bernoulli numbers Bj(0) = (— l) j Bj(l). 

Proposition 1. As n — > 00, we have C n (a + 1/2) — > —C n (a). 

Proposition 2. For integers m > and n > 1 we have 

7n = (-l) w -W-"X: 3(n + 1> ? +1 ~ fc) r P n {mx)^dx-J:C n {r/m). 

k=0 K - Jl X r=l 

Proposition 3. For integers g > 2 we have 
Proposition 4. For cr = Re s < 1 we have 

Proposition 5. For large but fixed fc, < a < 1, we have 

7fc(a) ~ msin[27r(a + 0)], 
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for real m and < < 1. 

Proposition 6. For integers n > 1 and k > 1 we have 



£C fc a± 



fc+i 



s(k + l,k + l - 



P k [n(x - a)]^r^rdx. 



x 



Proposition 7. For integers j > 1 we have the set of coupled differential equations 



(-l) 3 'rf7 3 -(a) 
j\ da 



k=j-i 



-l) fe /Jfe + l 



and 



-ip\a) = -l~J2 T\ 7fc(°)) 
fc=0 ft! 



. 7fe(o), j > 1, 



■1)* 



where ■?/> is the digamma function. 
Proposition 8. We have (a) 



1 1 ~ 1 " . k fn\\n(k+l) 



for the Euler constant, (b) 



-7i 



(c) 



»7i 



(d) 



In 2 2 1 ~ 1 



12 2 ^ 2 n+1 . , 

n=U fc=l 



'ra\ ln(k + 1) 



'n\ In 2 (A; + 1) 



V + 21n2 n ^ 2«+i^ V W + ' 



In 2 2 1 1 



12 + 21n 2 2 



h^ n+1 h K ) W + J ln2^ 2«+^i J W + ' 



ln7r = ln2-2£ — £(-l)*r In^ + l), 



On+l 
n=0 ^ fc=l 
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(e) 

(_1)Tt (2n)\ B2n = l-2 1 - 2 -§^ TT £ (_1)fc W + 

= (2^1) g JTI £ M)fe (*) + W n ^ °' 

and (f) the summation identity 

(^s^s ( - i K3 M4+i)( * +i)2 " 

^ ' 2(2ir) 2 "(l-2- 2 »)^- 2»+ 1 l t^ Wf^ + i) + 
In part (c), the constants r/j are denned by the expansion [7, 8, 9] ln£(s) — — ln(s — 

i) - e p °° =1 ^(s - iy, s ? i. 

These theorems are representative, but not comprehensive, of our results. We do 
not reproduce relevant background material on properties and functional equations of 
the special functions involved, but simply mention the useful sources [1, 13, 14, 26, 30]. 

Proofs of the Propositions 

Here we summarize or otherwise simply indicate the method of proof of the Propo- 
sitions above. 

Proposition I. A key observation is that 



_i2(2n)! 



P 2n ( x - a ) = (-l)"- 1 ^- [cos2tt(x - a) + 0(2~ 2n )\ , (1) 

and similarly for P2 n+ i(x — a), based upon the Fourier expansions of these polynomials 
[1]. So for large n, P2 n { x — a) = —P 2n (x — a — 1/2) + 0(2~ 2 ™) and we may appeal to 
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the representation [30] 

C n {a) = (-If" >n! J2 



11+1 s(n + l,n+ 1 - k) r°° „ , Jn fc x 



fc=0 



fc! 



/ p^x-a)-^^, n>l. (2) 

J 1 3/ 



Due to the boundedness of P„, \P n (x)\ < [3 + (—1)"]/ (27r) n for n > 1 [3], there is uni- 
form convergence of the integral in this equation. Therefore, we may use expressions 
like (1) in Eq. (2), interchange the limit and integration, and Proposition 1 follows. 
We have thereby proved a recent Conjecture II put forth by Kreminski [21]. 

Proposition 2. The multiplication formula satisfied by the Bernoulli polynomials 
[1] carries over to P n . Applying this to the representation (2) yields the Proposition. 

Proposition 3. One procedure is to take derivatives of the special case of the 
Hurwitz zeta function Y^1=i C( s ? r /<?) — (9 s — 1)C( S ) [14] and to evaluate as s — > 1 + . 
Another way is to use the defining Laurent expansions for the Stieltjes constants 
[21, 27, 30] in this special case, writing 

ge ( S -l)ln <? _ 1 l c(s) 



E 

r=l 



-l h k\ [ } 



k=0 



7=0 J ' 



k=0 



k\ 



(3) 



We expand the right side of this equation and cancel the term (q — l)/(s — 1) on both 
sides. We then reorder the double sum on the right side and equate coefficients of 
like powers of s — 1 on each side of the resulting equation. The Proposition follows 
once again. Proposition 3 takes the same form as Corollary 11 of Dilcher [10], who 
introduced a generalized digamma function. We explain this more in the Discussion 
section. 



Proposition 4. The key here is to note that 



f ((s,a)da = 0, 
Jo 



a < 1, 



(4) 



and then to integrate the Laurent expansion of the Hurwitz zeta function. Remark. 
Many other related integral-based results for the Stieltjes constants are possible. 

Proposition 5. For proof, we give the details for C 2n _i(a), with the case for C 2n (a) 
being very similar. From Eq. (2) we have 

\n k x 



C 2n -i{a) = (2n - 1)! 2^ - 

fc=0 



k\ 



J P 2n -i(x - a)—^-dx, n > 1. (5) 



We have [1] 



P (r n)-( ir 2(2n ~ 1)! 
F 2 n-i{x - a) - (-1) -j^-p^r 



sin27r(x-a) + 0(2 1 - 2n ) 



(6) 



so that we obtain 

^s(2n,2n-A;) 

^2n-lW — 2 / ~ 
fc=0 



fc! 



/ cos 2™ sin 27nr — sin 27ra cos 27r:r + 0(2 1_2n ) — - — dx, 
Ji L v y J x 2n 



(7) 



where «;„ = (— l) n 2[(2n— l)!] 2 / (27r) 2n x . As a result of the integration in this equation, 
the leading term in C 2ra _i(a) may be written as 



C2n-i( a ) ~ cos 2™ — r 2 (n) sin27ra, 



(8) 



where 



2n 



ri(n) = n n J2 



s(2n,2n - k) f°° ln k x 



k=0 



f°° in x , 

/ sm27rx — - — dx, 
Ji x 2n 



(9a) 



and 



r 2 (n) = K n Y^ 



k=0 



2n 



s(2n, 2n - k) 




(96) 



The form (9) can always be written as given in the Proposition with tan2-7T0 = —r 1 /r 2 
and m 2 = r\ + r\. Remark. The integrals appearing in Eqs. (9) can be written in 
several different ways via integration by parts, but we do not pursue this here. 

Proposition 6 follows from Eq. (2), the multiplication formula satisfied by P n [1], 
and the interchange of two finite sums. 

Proposition 7. According to the property d((s,a)/da = —s((s + l,a) and the 
Laurent expansion of £(s, a) we have 



wherein djo{d)/da = —ip\a) and ip' is the trigamma function. If we perform a 
binomial expansion on the right side of this equation, reorder the sums there, and 
then equate coefficients of like powers of s — 1 on both sides, we arrive at the stated 
set of equations. 

Proposition 8. Valid in the whole complex plane is the form of the Riemann zeta 
function 



This expression, due to Hasse [15], can be derived by applying Euler's series trans- 
formation to the alternating zeta function [25]. 




(10) 




(11) 
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For purposes of expanding Eq. (11) about s = 1 we have 

\n2) j 



3=1 J ' 



<s-iy, 



(k + ir s = T,^r ^(k + 1) E ( J (« - 



£=0 



(12a) 
(126) 



9=0 



and 



By using the series YZ=o (fc) fefi) = ^ * s easy ^° see ^ nat ^ ne 1 — term of 
Eq. (12b) contributes a ln2 term in Eq. (11): 



(13) 



Therefore, comparing with the Laurent expansion of £(s) we obtain parts (a) and (b) 
Similarly, if we write 

1 



lnC(s) = -ln(l-2 1 - s )+ln 



oo i n 

E^rE(-i) fc l 

n=0 ^ fc=0 



k (k + iy 



In 



= " H* - 1) + ^ In 2( S - 1) + 1( S - I) 2 + 0[( a - I) 4 ] 

E E (-D fc Q E ^ + 1) E Q (- - D< 

n=0 Z fc=0 W £=0 9=1 W 



(14) 



we find —i] = 7 and the expression given in part (c). 
From Eq. (11) we have 

In 2 1 ~ 1 ™ ,Jn\\n(k+ 1) , , . 
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Part (d) follows from the well known value ('(0)/((0) = \u2tt. Similarly, based upon 
the well known values of the zeta function at positive even integers we obtain part 
(e) for the Bernoulli numbers. 

From the functional equation of the zeta function comes the evaluation for positive 
integers n 

«-*»> = w 

From the representation (11) and another globally but more slowly convergent series 
for ((s) due to Hasse [15] we then have the identity of part (f). 

Discussion 

We may note the special values of the Stirling numbers of the first kind [1] s(n, 0) = 
5 n o and s(n, n) = 1. Therefore, the general integral representation of Eq. (2) can be 
simplified to 

P n (x - a) ™ s(n + + 1 - k) f°° ln k x 



C n (a) = {-l) n - l n\ 



r^^ n {x-a "s{n + L,n + L-k f°° 



, dx 

, , k\ Ji ' x n+1 

(17) 



When a — 1, the first term on the right side of this equation may be evaluated by 
using integration by parts on the result of the Appendix. 

Despite some of the complexity of the representation (11), we have several motiva- 
tions for examining it. (i) This form gives rise to series that play a significant role in 
renormalization and quantum field theory [5, 12]. (ii) Equation (11) has connections 
to random variables in an analytic number theory setting [4]. There is a possibility 
in particular to link this structure with that of Brownian processes. This brings us 
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n > 1. 



to a third point of interest, (iii) We would like to know not only the structure of the 
Stieltjes constants but also of the constants {rjj} that enter the Laurent expansion of 
the logarithmic derivative of ((s). Such knowledge would be very influential in decid- 
ing the nature of a sum denoted S 2 (n) [7, 8, 9] and thereby the resulting Li/Keiper 
constants [22, 19]. 

As a byproduct of this work we obtain interesting infinite series for fundamental 
constants such as the Euler constant and ln7r. The rapidity of convergence may make 
some of these suitable for applications. Indeed, even a naive numerical implemen- 
tation of part (a) or (b) of Proposition 8 appears to provide these constants to 16 
decimal places after summing over n to 51 or 52 terms. 

The generalized digamma function used by Dilcher [10] is given by 



■y OO 

ipk(a) = -Ik hi fc a - 



lnV + a)_ln^l 
v + a v 

and we may relate it to other functions defined in the literature. We have R m (a) = 
(-l) m+1 (d m /ds m )((0,a) [17] and find that ip k (a) = R' k+1 (a)/(k+l). We have Rj(a) = 
{—\yj\ — J2tLo lj+i( a )/^- an d by using our Proposition 7 we determine that ipj(a) = 
— 7j(a). This relation explains how our Proposition 3 takes exactly the same form as 
Dilcher 's Corollary 11 and permits the re-expression of several of his other results as 
well. 

In his computational work based upon Newton-Cotes integration for the high 
accuracy approximation of the Stieltjes constants, Kreminski observed that for large 
values of k, 7 fc (l/2) ps — ^ k and that more generally C k (a + 1/2) ps —C k (a). That is, 
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in a sense, the Stieltjes constants for large index are anti-periodic with period 1/2. 
Our Proposition 1 makes this precise. 

Much earlier, Hansen and Patrick [14] showed that a fundamental interval in a for 
the Hurwitz zeta function need only be of length 1/2. Perhaps solely on that basis 
one would then suspect that some sort of relationship (s) should exist between 7^(0) 
and 7 fe (a + 1/2). 

Much of our development relies on the underlying theory of periodized Bernoulli 
polynomials and corresponding integral representations of the Stieltjes constants [30]. 
Equations such as (17) and still others that we have obtained help to expose more of 
the analytic structure of the Stieltjes and rjj constants. 
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Appendix: Integrals over periodic Bernoulli polynomials 

Herein we evaluate integrals over the periodic Bernoulli polynomial P\ in terms 
of polygamma functions ip^\ We demonstrate 
Proposition for positive integers n and m 



-l) n 



(n + 1) 



1 



E 



i 



(n + 1) f & + k) n+1 
1 1 



Pi(x) 



-dx+ 



1 



1 



x + s + l) n + 2 (n + 1) (s + m) n+1 

11 1 /A \ 

(Al) 



n(n + 1) (s + m + l) n 2 (n + 1) (s + m + l) n+1 
Our starting point is [11] 

Pi(x) 



i r°° 
lnr(s + l) = (s + l/2)lns-s + -ln27r 

2 



Jo 



(x + s) 



dx, 



so that 



1 f c 
^(s + 1) = ms + — + / 
is Jo 



Pi(x) 



dx, 



(A2) 



(A3) 



(x + s) 2 

where ip(s) = ip(s + 1) — 1/s is the digamma function. We now take n derivatives of 
relation (A. 3), obtaining 



(-1)" 



(n + 1)! 
Finding that 
Pi(x) 



^ (n) (s+l) 



1 1 1 

+ 



+ 



n(n + l)s n 2 (n + 1) s n+1 Jo (x + s) n+2 



Pi(x) 



-dx. (A4) 



-dx = 



dx+ 



o (x + s) n + 2 ' Ji (x + s) n + 2 n(n + 1) 



(s + l) n 



2(n + l) 



(s + l) n 



+1 



(A5) 



we obtain 

(-1)" 
(n + 1) 



1 1 



(n+l)s n+1 n{n + l){s + l) n 2 (n + 1) (s + 1) 



l r c 

i) n + 1+ Ji 



Pi(x) 



(x + s) n + 2 
(A6) 
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By making a simple change of variable in the integral of Eq. (A. 6) and using the 
periodicity of Pi, we have 

D" (B) , 1 1 1 1 



^ n) (s + m + 1) 



(n + 1)! ' ' (n+ 1) (s + m) n+1 n(ra + 1) (s + m + l) n 

1 -fr^U- <-) 



2 (n + 1) (s + m + im (x + s + l) n+2 

Applying the functional equation of the polygamma function [1] we obtain the Propo- 
sition. 

When applying Eq. (A.l), it is useful to keep in mind a relation [13, 26] between 
the polygamma function and the Hurwitz zeta function: ^ n \x) = (— l) n+1 n!£(n + 
l,x). In particular, when x = 1/2 or x = 1 in the latter relation, values of the 
Riemann zeta function appear. In fact, another starting point for evaluating integrals 
similar to (A.l) is to use integral representations for ((s) [1, 16, 26, 28] in terms 
P 1 (x) + l/2. 
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